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Abstract 

A theory describing above-threshold ionization of atoms and ions in a 
strong electromagnetic field is presented. It is based on the widely known 
strong field approximation and incorporates the Coulomb interaction be- 
tween the photoelectron and the nucleus using the method of complex clas- 
sical trajectories. A central result of the theory is the Coulomb-corrected 
ionization amplitude whose evaluation requires little extra numerical ef- 
fort. By comparing our predictions with the results of ab initio numerical 
solutions for two examples we show that the new theory provides a signif- 
icant improvement of the Coulomb-free strong field approximation. For 
the case of above-threshold ionization in elliptically polarized fields a com- 
parison with available experimental data is also presented. 



1 Introduction 

Among analytical approaches developed for the description of strong field ion- 
ization of atoms, ions and molecules, the so-called strong field approximation 
(SFA) or Keldysh-Faisal-Reiss theory (KFR) [TJ [U [3] is known to be the most 
fruitful and widely applicable. The objective of SFA is the description of above- 
threshold ionization (ATI) of atoms and ions in a strong electromagnetic field. 
Being generalized to account for the effects of rescattering, SFA was used to cal- 
culate the high-energy plateau in ATI spectra, high-order harmonic generation 
and nonsequential multiple ionization. The present state of the theory is re- 
viewed in Refs. $, 5J El U\ ■ During the last years the scope of strong field physics 
has been extended to more complex systems including molecules, fullerenes and 
clusters where the SFA is also being used as a powerful tool (for SFA in molecules 
see, e.g., Ref.[5] and references therein). 

There are several formulations of the SFA which are not fully equivalent. 
However, all of them are based on the idea that the electron continuum states 
can be approximated by plane Volkov waves. The latter are exact solutions of 
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the Schrodinger equation for a free electron in the field of a plane electromag- 
netic wave [5] . Using the plane Volkov continuum the matrix elements describing 
ionization and other related processes admit fast numerical evaluation. More- 
over, simple analytical treatments are possible in limiting cases. Physically, the 
SFA fully disregards the effect of the binding potential on the continuum. For 
negative ions (or, more general, for systems where electrons are bound by short- 
range forces) , where the Coulomb interaction between the detached electron and 
the atomic core is absent, this approximation is well justified. As a consequence, 
the SFA usually provides a quantitatively correct description of ATI and related 
processes in negative ions, as has been multiply proved by comparing SFA re- 
sults with data [TOl HU H2] , with predictions of a more rigorous quasistationary 
quasienergy states method [T31 Q31 HH [H] and with numerical solutions of the 
time-dependent Schrodinger equation (TDSE) for single active electron model 
atoms [T7] . 

In case of atoms, positive ions and molecules the neglect of the long-range 
Coulomb interaction between the emitted electron and the ion core is not justi- 
fied. On a qualitative level the SFA often works for such systems too. Several 
SFA calculations also led to good quantitative agreement with experimental data 
[18] . However, it is now well established that the neglect of the Coulomb inter- 
action may lead to essential contradictions. Violation of the gauge invariance 
and wrong predictions for the static-field limit in the total (energy-angular in- 
tegrated) ionization rate are widely known examples. Besides, there are several 
well-documented features in ATI spectra which certainly cannot be explained 
within the standard SFA approach (see [HSE!] for examples). 

In recent years new opportunities to clarify the quality of different theoretical 
approaches are possible by ab initio numerical solutions of the time-dependent 
Schrodinger equation for one and two-electron systems driven by strong laser 
fields. Compared to real experiments these numerical experiments provide more 
stringent tests for analytical theories since all parameters are exactly known and 
no focal averaging washes out interference effects which are sensitive to the de- 
tails of the laser pulse shape and the carrier-envelope phase [21] . In general, 
the results of numerical experiments show that the SFA is only accurate on 
a qualitative level in the description of ATI angular-resolved spectra, even in 
the simplest case of atomic hydrogen. Many important details, including the 
interference structure, symmetry properties of photoelectron momentum distri- 
butions and the positions of ATI maxima reveal a disagreement between the 
SFA and TDSE results |T7l |22] . The key physical effect behind these disagree- 
ments is the long-range Coulomb interaction between the emitted electron and 
the atomic core. 

Since the early days of laser physics considerable effort has been invested into 
taking the Coulomb interaction in strong field theories into account. The first 
attempt was probably made in 1967 by Nikishov and Ritus |23j . In the same 
year, a Coulomb correction to the tunnel ionization rate of atoms was derived by 
Perelomov and Popov [21] . In the 80s the Coulomb- induced orders-in- magnitude 
enhancement of tunnel ionization rates of atoms and positive ions was well doc- 
umented in experiments [25 . Later, different models were proposed to account 
for the effect of the Coulomb force on the photoelectron during its propaga- 
tion from the atom to a detector within classical mechanics. Such calculations 
explained to some extend the Coulomb-induced asymmetries in elliptically |26j 
and linearly [27j polarized fields; in the latter case an asymmetry may only ap- 



2 



pear in very short pulses [2T]. Classical modeling of photoelectron trajectories 
in the presence of both the laser and the Coulomb fields has shown peculiarities 
in the extremely-low-energy part of ATI spectra H§] . These peculiarities 
in the form of cusps and dips in photoelectron momentum distributions were 
recently observed in high-resolution experiments [30\ 120] . Within quantum ap- 
proaches, so-called Coulomb- Volkov functions have been used for the description 
of Coulomb effects in ATI [3TJ [33J (Ml HHI- For molecules, these wave functions 
were generalized in [31] . In Ref. [35] the Coulomb-phase-distorted wave function 
was proposed and applied to the analysis of momentum distributions recorded 
in [30l [20] . Finally, in Refs.[36l [37] a new eikonal-like approximation for the 
Coulomb- Volkov continuum was formulated and applied to the description of 
the subcycle ionization dynamics and the strong-field assisted XUV ionization. 
However, to the best of our knowledge up to now no reasonably simple the- 
ory of strong field ionization incorporating the Coulomb field in all its essential 
features has been put forward. 

In this paper, we show how the SFA should be modified to incorporate 
Coulomb effects. We call this new approach the Coulomb-corrected strong field 
approximation (CCSFA). The basic idea is to account for the Coulomb field 
by using semiclassical perturbation theory for the action. In connection with 
the problem we consider this method was proposed in 1967 by Perelomov and 
Popov [24] . Here we generalize this idea for the description of angular-resolved 
ATI spectra. This generalization exploits the imaginary time method [35] or, 
equivalently, the technique of complex trajectories [321 SOI HI] • Applications of 
the proposed method demonstrate significant improvements of the theory. In 
fact, its predictions appear to be almost identical to the corresponding results 
of the ab initio single electron TDSE solution and yields quantitative agreement 
with experimental data. Although the theory requires a few more computations 
than the standard SFA, it remains remarkably simple. 

The paper is organized as follows. The next section is devoted to the stan- 
dard SFA ionization amplitude and to the semiclassical picture of complex clas- 
sical trajectories. In the third section a general method to include Coulomb 
corrections is described and the CCSFA ionization amplitude is formulated. 
The subsequent section shows two applications of the new theory. The last sec- 
tion contains final remarks and conclusions. Atomic units are used unless noted 
otherwise. 

2 Coulomb-free Strong Field Approximation 

2.1 SFA matrix element 

Within the SFA the transition amplitude between an atomic bound state |^o) of 
binding energy eo = —I and a continuum state \^ p ) with asymptotic momentum 
p is given by 




(1) 



— OO 
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where the final state is replaced by the Volkov function 



* p (r, t) = exp j ipv - i ^ (p + A^fdh I (2) 



and V"(i) is the interaction operator which may be chosen in different gauges. 
The laser field is described by the vector potential A(t) which depends only 
on time in the dipole approximation we use here. The electric field is given by 
£{t) = -d t A. 

For a ground state in a short-range potential the amplitude ([T|) is gauge- 
independent, so that we equally may use the Volkov functions and the inter- 
action operators in the length or velocity gauge. If one starts from an atomic 
bound state Msfa(p) in general is gauge-dependent, with sometimes very sig- 
nificant differences between length and velocity gauge. For atoms the length 
gauge was shown to be the better choice [HI [12] while for spatially extended 
systems such as molecules the velocity gauge works better [§]■ However, in 
our following considerations we shall not be confronted with this unphysical 
gauge-dependence. We start from the matrix element |T]) for a state in a short- 
range well of ionization potential I equal to the one of the real atomic state we 
consider. 

Two dimensionless parameters (defined differently in the literature) turn out 
to be important within the SFA. Here we use the Keldysh parameter 7 and the 
multiquantum parameter K$, defined as 

7 = kuj/£ , K = I/uj (3) 

where £q and u> are the field amplitude and the laser frequency, respectively, 
and k = V2I is the characteristic momentum of the bound state. 

The parameter Kq defines the minimum number of photons necessary for 
ionization. If Kq 3> 1 then, at arbitrary 7, the integral over time in J!} may be 
evaluated by the saddle-point method so that the amplitude can be represented 
as a sum of contributions from all relevant stationary points t s (p) , 

M ( n \ ri)/nf , exp{-iS (p,t sa )) 
M SFA (p)-X>(p,W 7Wma » W 

where 

+00 

S (p,t)= j |i(p + A(ti)) a + l|dti (5) 
t 

and the pre-exponential V contains the spatial matrix element. The saddle-point 
equation has the form 

d t S (p,t sa ) = i(p + A(t SQ )) 2 + / = (6) 

which shows that a saddle point t s (p) is always complex for I > 0. The differ- 
ential ionization rate is given by the square modulus of J!]). 
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2.2 Formulation in terms of complex trajectories 

The amplitude Q can be equivalently reformulated in terms of classical complex 
trajectories [SB]. For each stationary point consider a trajectory which satisfies 
the Newton equation in the laser field (below we omit a subscript a unless it is 
misleading) 

r = v = -£{t) (7) 
with the initial and boundary conditions 



v£(i = t a ) = -k\ v (i^+oo) = p. (8) 



The solution is 



t 



r (p,i)=Xo + p(i-ia) + G(t)-G(t a ), G(t) = J A{ti)dt\. (9) 

Although a trajectory satisfies the classical equation of motion it is in general 
complex because of the initial condition. However, the velocity vo = p + A(t) 
is always real in real time. We assume that the real part of a trajectory starts 
from the position of the atom, Re[ro(p, t s )] = 0. Then, by a proper choice of 
a purely imaginary constant Xo, we can make the coordinate real in real time 
too. The trajectory ^ starts from Xo at t = t s , having a purely imaginary 
velocity, Vg(p,£ s ) = — k 2 . As time approaches the real axis, t = to = Rc[i s ], 
both velocity and coordinate become real, so that the electron appears in real 
space at 

Ro(p) = r (p,to) = G(i )-Re[G(t s )]. (10) 

The position Ro is commonly interpreted as the "tunnel exit". Although the 
tunneling picture of ionization is only relevant for 7 <C 1, the trajectory- 
approach is valid at arbitrary values of the Keldysh parameter. Hence, we 
will use the term "exit" for (jTDJ) for arbitrary 7. Beyond the exit, the elec- 
tron propagates in real space and time towards a detector. The values Ro and 
Vo = Vo(p,to) should be considered as initial conditions for the motion in real 
time. They fully determine the real part of a trajectory. 
After the following algebra, 

+00 

So(p)= J Qvo+/+|(vo-ro)-|(vo.r ))d^ 

t.(p) 

= - J Qvq - v ■ r - l\ dt + (v ■ r )|t^+ oo - (v ■ ro)| t = t ,, 

ts(p) 

we represent (j4]) in the form 

M SFA (p)=E^(p) f ^p# (id 

where 

+00 

W (p)= J (£ -I)dt-p-r Q (+oo)+v {t s )-X (12) 
*.(p) 
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is the classical reduced action evaluated along the complex trajectory @. The 
corresponding Lagrangian is Co — Vq/2 — £ (t) • r . 

The result (fTTj) shows that the SFA ionization amplitude can be represented 
as a coherent sum of contributions from complex classical trajectories. To cal- 
culate each contribution we have to find a trajectory satisfying Eqs.((7|), © and 
evaluate the respective classical action (TT2"j) . 



3 Coulomb-corrected Strong Field Approxima- 
tion 

3.1 Basic idea 

The above formulated representation of the SFA provides a natural way to in- 
clude the Coulomb field via corrections to the classical action in [fTTj) evaluated 
along classical complex trajectories. Assuming these corrections to be small 
compared to the laser-induced action (|12jl we may perform this procedure per- 
turbatively taking Coulomb-free trajectories as a zero-order approximation [24j . 
There are two corrections. One is due to the appearance of the potential energy 
Uc in the action 



+oc +00 



wg 



(P) = -/ U c[ r oM]d t = Z J * , (13) 



< s (p) *»(p) 
another one is due to the correction ri to a trajectory 

+00 

W^ U \p) = J [v • V! - £(t) ■ n]di. (14) 

Here the correction ri(p, t) has to be determined from the Newton equation 

2(r + ri) 

1-1 = — i ! iF ( 15 

with Z being the residual charge (Z = 1 for a neutral atom). Besides, the 
stationary point t s (p) may also be corrected. 

This general idea meets, however, with several difficulties: 

(i) Corrections ([13)) and (Til)) are logarithmically divergent at the lower integra- 
tion limit t — > t s when the electron approaches the nucleus because for the most 
important trajectories the imaginary starting point Xo in ^ is small or exactly 
equal to zero. Obviously, under such circumstances the singular Coulomb force 
cannot be treated perturbatively. 

(ii) A zero-order trajectory i"o(p, t) may revisit the nucleus in real time (such 
trajectories are known to be responsible for the rescattering phenomena). This 
causes another divergence. 

(iii) It is unclear which initial conditions should be attached to Eq. (fT5|) for the 
correction ri(p, t). 

The problem (i) was solved already in Ref. [21] • Below we show how to resolve 
the other two by the proper choice of a zero-order trajectory. 
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3.2 A choice of the zero-order trajectory 

If we take a trajectory ro(p,t) and evaluate a correction in real time by solv- 
ing (|T5|) for t > to then, even if the electron does not revisit the origin and the 
correction is everywhere small, we obtain a different final photoelectron momen- 
tum v(t — ► +00) 7^ p. Since, however, we want to find the ionization amplitude 
for a given final momentum p, the above described procedure does not appear 
satisfactory. 

A correct realization is as follows. Assume that we know some other Coulomb- 
free trajectory i"o(p, t) such that, if we account for the Coulomb force at t > to, 
the final momentum will exactly have the desired value p. Below we denote 
such a trajectory r (p,<) = r (p,t). A corresponding stationary point t s (p) 
satisfies Eq.© with p instead of p so that 

?.(p) = i.(p). (16) 

The new initial conditions for the motion in real time have the form 

R (p) = R (p) = G(i ) - ReG(t.), V„(p) = p + A(t ). (17) 

How to find this trajectory? First, we give some physically transparent ar- 
guments which explain why such a trajectory and the respective renormalized 
momentum p always exist. We start the electron's propagation at t — to when 
the electron is at the "tunnel exit" Ro- Here the Coulomb force is already small 
compared to the laser field (for an exact criterion see Subsection 3.4). Therefore, 
if the electron does not approach the nucleus, a Coulomb-induced deformation 
of the trajectory is smooth. Thus, just according to physical continuity each tra- 
jectory ro(p, t) must have an analog ?o(p, t) which, without the Coulomb force 
corresponds to some different final momentum p and goes to p if the Coulomb 
field is turned on. Moreover, even if the initial trajectory is approaching the 
nucleus, the new one will not because of the following two reasons. 

1. In full dimensionality the fraction of trajectories which closely encounter the 
nucleus is very small. Even in a linearly polarized field the relative amount of 
trajectories which approach the parent ion at distances of several atomic units 
typically is 1CP 4 -j- 10~ 6 . This value is nothing else but the relative probability 
of rescattering, known to be always small. 

2. If a trajectory approaches the nucleus, it appears to be strongly disturbed 
by the Coulomb force and, as a result, the electron in general will contribute to 
the high-energy rescattering part of the spectrum we are not interested in here. 
Therefore, such trajectories, although existing, cannot be those we are looking 
for. 

In a more strict formulation, although some fraction of trajectories appear 
to be strongly disturbed by the Coulomb field and show irregular behavior, 
for any given value of the final momentum p new trajectories fo(p, i) exist 
which, after the Coulomb field is taken into account, experience only a smooth 
perturbation. The number of such regular trajectories is equal to the number 
of initial trajectories ro(p, t). 

We illustrate this statement for the case of a linearly polarized field where 
the most important trajectories may approach the nucleus closely so that this 
example is the most challenging one. Consider ionization in the field 

£(t)=£ cos(ut), (18) 
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x, a.u. x, a.u. 

Figure 1: (Color online) Panel (a): Trajectories corresponding to ionization from 
the ground state of hydrogen into the final state with energy e = p 2 /2 = 0.278 
(5-th ATI maximum) and emission angle 9 = 0.05 in the field (fT5|) with £q = 0.08 
and lo = 0.056. The emission angle is chosen different from zero just to make the 
Coulomb- free trajectories ro±(p, t) two-dimensional and therefore visible in the 
plot. They are shown by green and black solid lines, the respective trajectories 
fo±(p) by dashed lines. The full solutions of Eq. (fT5|) r±(p, t) with initial condi- 
tions (fT7|) are shown by red and blue. An irregular trajectory developing from 
r - under the action of the Coulomb force is shown in gray. Panel (b) shows the 
space indicated by a red rectangle in panel (a) with more resolution. The posi- 
tion of the atom is indicated by a red circle. It is seen that the trajectory which 
does not return experiences only a small distortion while the other trajectory 
is changed significantly. The values of the final and the renormalized momenta 
are p = (p x ,p y ) = (0.745,0.037), p+ = (0.971,0.043), p = (0.369,0.223). 
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polarized along the x-axis. For a constant field envelope all laser periods are 
identical. For each final momentum p there are two trajectories per period (we 
denote them i"o+(p, t) and ro_(p,i)), whose motion in real time starts from 
the points Ro+ and Ro- = — Ro+ symmetrically displaced with respect to the 
atom. If the transverse momentum p y is small or equal to zero the trajectory 
ro+ will never return while ro- revisits the nucleus at least once (see Fig.f). It 
is intuitively clear that the Coulomb force along the first trajectory is always 
small and simply decelerates the electron on its way to a detector. Therefore, 
the new Coulomb-free trajectory fo+ will be quite close to the initial one and 
the renormalized momentum p + is expected to be larger in absolute value than 
p in order to compensate this deceleration. In contrary, the electron moving 
along ro_ will accumulate a significant negative transverse momentum due to 
the Coulomb attraction when it passes near the nucleus. Using ro_ as a zero 
approximation we obtain a typical irregular trajectory shown in Fig.l in gray. 
Thus, to arrive at a detector with a prescribed small transverse momentum, 
the electron has to start with some significant one, p y > 0. However this 
immediately leads to the trajectory fo_ which stays away from the nucleus 
and remains regular. This example, summarized in Fig.l, shows that even in a 
linearly polarized field the regular Coulomb corrected trajectories exist and they 
can be obtained from initial Coulomb-free trajectories by a smooth deformation. 
The key point is to find proper zero-order trajectories fo instead of ro. In an 
clliptically polarized field, where trajectories are always two-dimensional, the 
situation is even simpler. 

Another question is how to find these trajectories in a particular calculation. 
Except limiting cases this can only be done by solving Eq. IT5|) numerically. Since 
not the initial but the final condition v(< — > +00) = p is fixed we have to deal 
with an inverse problem. Because of this additional complication a proper tra- 
jectory can only be found, to the best of our knowledge, by iterations starting 
from some appropriate zero-order approximation. However, this numerical it- 
eration is simple to implement and the evaluation of one Coulomb-corrected 
trajectory takes a fraction of a second on a PC using a standard solver for dif- 
ferential equations. In the two applications shown below two different iteration 
procedures have been applied, both leading rapidly to converged results. 

3.3 Sub-barrier corrections and the matching procedure 

The algorithm described in the previous subsection is incomplete since correc- 
tions during the "sub-barrier" motion while i 6 [t s , to] should also be taken 
into account. Because a trajectory fo already yields, after accounting for the 
Coulomb force at t > t s , the desired final momentum p, the boundary con- 
ditions for sub-barrier motion are ri(tg) = 0, Vi(to) = 0. Consider the most 
"dangerous" case when Xo = and all the values ri, vi, W^p and Wq are 
divergent at t — ► t s . The regularization procedure is based on the idea that, as 
the electron approaches the nucleus, its dynamics should be determined mainly 



by the binding force, not by the laser field [24 . Thus, the respective action W&t 




can be found from an asymptotic expression for the atomic bound state wave 
function at large distances r S> 1/k. For s-states this asymptotic behavior is 

^o( r ) ~ («t)™* exp{— «?•} = exp{— ft?' + n* ln(ftr)} = enp{iW a t(r)} , 
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so that 

M / at (r) = i[nr — n* ln(rer)]. (19) 
Here, n* = -E/k is the effective principle quantum number of the bound state 

To match the Coulomb-corrected action with (fP9|) we extract the divergent 
part of the former at t — > t s . To this end we introduce a matching point 
t* £ [t s ,t ] such that 1/k < |r (p,**) +ri(p,i*) = r*| < i?o(p)- Here, the 
laser-induced action (TT2"]) has the behavior 

W (t*)« W (p) - 2/(t. - f). (20) 

Integrating (JT4j) by parts and taking into account ([7]) and the initial conditions 
r i(^o) = v i(^o) = we obtain 

W$ I) (t*) = -iK(e(t s )-r* 1 ), e(t s )=£(t s )/\£(t s )\. (21) 

Finally, the divergent part of (jT^J) is 

Wg\t*) m -in. In (22) 

At i* — > t s the solution of Eq.© with the initial conditions ro(t s ) = 0, Vg(i s ) = 
— k 2 is 

r*Er (t>-ke(i 8 )(t,-f). (23) 

Now, collecting the expressions ((HI), PD}. {3D]), and fl2TJ and pg]). we see that 
indeed all terms depending on the matching point disappear in the full action 
Wat(r*) + W (t*) + W§\t*) + W£ l \t*) and the final result is convergent. 

The regular part of p^|) is equal to zero, the one of (fTB")) can be evaluated for 
arbitrary field polarization and photoelectron final momentum only numerically. 
In a linearly polarized field, for the most probable trajectory with p = 0, it is 

m 

Reg[VF C i](p = 0) = -m« ln(2K 3 /£o) = -in* ln(2/F), 

where F — £o/k 3 is the reduced laser field; the value £ at = k 3 has the meaning of 
a characteristic atomic field. This contribution causes the well-known Coulomb 
correction enhancing the total ionization rate [23] . The regularization procedure 
described above becomes more cumbersome, although doable, for a laser field of 
arbitrary polarization and for an arbitrary final momentum if Xo(p) ^ 0. Also 
bound states of arbitrary angular momentum can be treated. 



3.4 CCSFA ionization amplitude and applicability condi- 
tions 

Summarizing the results obtained in this section we formulate an algorithm for 
the evaluation of the CCSFA ionization amplitude: 

1. For a given value of the final momentum p determine all relevant saddle 
points i sa (p) and the respective trajectories r 0a (p,t). 

2. For each trajectory find the renormalized Coulomb- free trajectory ro a , the 
stationary point t aa , the momentum p a and the Coulomb corrected trajectory 
r a which is a solution of the Newton equation in two fields with the initial 
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conditions (JTTJ) by solving the inverse problem using an iteration scheme. 

3. Calculate the Coulomb- free action (fl^|) for the new trajectory fo a . Here 
one should note that the renormalized momentum p a enters only the imaginary 
part of the action since the latter is accumulated for t £ [t sa , toa], while the real 
part of the action depends on the final momentum p. 

4. Calculate the regular part of the Coulomb-induced action (TT3")) along the 
exact trajectory r Q . 

Finally, the CCSFA ionization amplitude has the form 

exp(i[W 0a (p) + W§l(p)}) 

where 

Im[Wo(p)] = Im[Wo(p,? s )], Re[Wo(p)] = Re[W (p,t s )]. 

Since the pre-exponential is much less significant we neglect corrections to it in 
this work. 

The amplitude is relevant only if the formalism of trajectories is appli- 
cable. This requires the validity of the saddle-point method for the evaluation 
of the Coulomb- free amplitude (P), i.e. Kq >■ 1. Moreover, the Coulomb force 
at the "exit" (JTUJ) must be small compared to the laser field amplitude. Using 
the reduced field F and the Keldysh parameter 7 the latter restriction can be 
formulated as 

1 > F x ( \ 7 < \ . (25) 

The condition (|25[) shows that the standard SFA can only be used as a zeroth- 
order approximation for the development of a quantitatively correct theory of 
strong field ionization in the case of relatively weak fields £q <C £ a t of not too 
high frequency u> -C IyfF. In terms of the laser field amplitude, the condition 
(|25|) can be reformulated as 1/Kq < F < 1. The corresponding intensity 
domain does not coincide with the one determining the tunneling regime, 7^1. 
The upper restriction on the Keldysh parameter 7 <C 1/ V~F is much softer and 
allows treatments in the multiphoton regime. Of course, for F — > or 7 — > 00 
condition will be violated, so that the amplitude (JH|) does not match the 
perturbation theory result. This feature is known to be inherent in the standard 
SFA as well. 



4 Applications 

In this section we evaluate the CCSFA ionization amplitude (|24[) and calculate 
the respective spectra for elliptically and linearly polarized fields and compare 
our results with the predictions of the standard SFA and those obtained by 
numerically solving the time-dependent Schrodingcr equation in three spatial 
dimensions for atomic hydrogen, neon and argon using the code described in 
|42| . The TDSE is solved for an atom with a single active electron whose bound 
state is characterized by the binding energy, orbital and magnetic quantum 
numbers. The self-consistent effective single electron potential U(r) may be 
constructed using the density functional method [43] or approximating it by a 
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screened Coulomb potential with parameters adjusted to reproduce both the 
correct asymptotic behavior U(r 3> 1/rt) ~ Uc{r) = — 1/r and the ionization 
potential of the ground state. The pulse with the carrier frequency u>, amplitude 
£o, ellipticity £ and a duration of n = 8 optical cycles polarized in the (x, y)- 
plane is described by the vector potential 



A(t) 



(-) 

\2nJ 



cos tp, 



COS if 



(26) 



with tp = u)t, tp' = tp — 7r/2 and < tp, tp' < 2irn. 



4.1 Symmetry violation in an elliptically polarized field 

It is now well established that, while the SFA predicts a fourfold symmetry of 
the photoelectron distribution for an arbitrarily polarized laser field (if finite 
pulse duration-effects are neglected), the experiments show that in elliptically 
polarized laser light this symmetry is clearly broken [19l 04]. According to 
the intuitive semiclassical ionization picture [45] , which often is in remarkable 
agreement with the predictions of the SFA, the most probable photoelectron 
momentum is along the minor polarization axis of the polarization ellipse for 
intermediate ellipticities, say, 0.2 < £ < 0.8. This effect of classical dodging 
[46] [47] simply follows from the fact that the electron tunnels out with zero 
velocity at the "tunnel exit" with highest probability near the maximum of the 
electric field. At that time the vector potential is oriented along the minor 
polarization axis, determining the drift photoelectron momentum p = v — A 
measured at a detector. The SFA angular distributions shown in Fig. 2 by green 
lines clearly demonstrate this effect. However, the experimental results do not 
display the fourfold symmetry at all: photoelectrons go predominantly along 
the major polarization axis and the distribution possesses inversion symmetry 
only. A typical distribution recorded in [44j is shown in Fig. 2(d). 

Already in early works [1^1 [H] it was suggested that the asymptotic Coulomb 
interaction of the outgoing photoelectron with its parent ion breaks the four- 
fold symmetry. However, except [26], no explicit theory- us-experiment analysis 
was presented. Although the results obtained in Ref.[55] demonstrate reason- 
able agreement with the data, the lack of ab initio calculations did not permit 
to draw definite conclusions. Now we may present the full analysis along the 
line standard SFA vs CCSFA vs ab initio TDSE results vs experiment. The 
iterative procedure in CCSFA was performed taking the trajectories ro(p, i) as 
a zero approximation and turning the nuclear charge on smoothly. Figure 2 
shows polar plots of the angular distributions for fixed ATI maxima calculated 
for different atoms and field parameters. Panel (d) shows the data for argon 
recorded in Ref. [44j . The standard SFA with its fourfold symmetry is in strong 
disagreement with both experimental data and numerical results. Agreement 
between CCSFA and ab initio TDSE results is everywhere quantitatively good, 
particularly for the cases (a,c) where the parameter ([23]) is smaller. The agree- 
ment with the data is also satisfactory, although neither CCSFA nor ab initio 
calculations (the latter is not shown in Fig. 2(d) because it is almost identical to 
the CCSFA result) show the well-developed second maximum which is clearly 
visible in the data. Among other possible reasons for this disagreement we think 
it is due to the uncertainty in the experimental laser intensity. The intensity de- 
clared in Ref.[H] was 0.9 • 10 14 W/cm 2 . In calculations we used 0.6 • 10 14 W/cm 2 
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Figure 2: (Color online) Normalized angular distributions in the polarization 
plane, evaluated by different methods: the standard SFA (green), CCSFA (red) 
and the ab initio TDSE solution (black) for ground states of hydrogen (a,b), 
neon (c) and argon (d) . In the latter case the data recorded in [H] are shown by 
black circles while the TDSE result is not shown. The laser intensity, frequency, 
ellipticity and the photoelectron energy are: 1.0 • 10 14 W/cm 2 , 1.55eV, 0.5 and 
8.2eV (a), 1.6 • 10 14 W/cm 2 , 3.1eV, 0.5 and 2.6eV (b), 2.0 • 10 14 W/cm 2 , 1.55eV, 
0.36 and 7.1eV (c) and 0.6 • 10 14 W/cm 2 , 1.55eV, 0.36 and 2.8eV (d). For argon 
and neon the ionization probability is averaged over the magnetic quantum 
numbers m = 0, ±1 within the p-shell. The orientation of the polarization 
ellipse is shown in the insert with the rotation direction of the electric field 
vector indicated by an arrow. 
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which led to better agreement including the signature of the second maximum 
(red curve). 

4.2 Interference in a linearly polarized field 

In a field of linear polarization, at least two trajectories per laser period con- 
tribute to the amplitude with comparable weights (exactly equal for the field 
(fTSjl having a constant envelope). A number of relevant trajectories > 1 leads to 
the interference structure visible in ATI spectra and was first considered in [35] . 
With increasing field cllipticity, for the vast majority of final momenta, the con- 
tributions from the above mentioned trajectories become significantly different 
in absolute value, and interferences disappear. In short pulses, the interference 
structure depends also upon the absolute phase, as it was observed in the at- 
tosecond double-slit experiment [49]. Ab initio TDSE calculations performed in 
Ref. [3H] showed good quantitative agreement with the data. The SFA results, 
however, reproduced the interference structure only on a qualitative level. An 
example shown in Fig. 3 demonstrates the obvious quantitative disagreement 
between TDSE and SFA: interference maxima and minima in SFA and TDSE 
results are almost opposite to each other, meaning that the interfering trajec- 
tories in plain SFA have a relative phase which is wrong by n. 

The spectrum calculated using the CCSFA (red line in Fig. 3) demonstrates 
good quantitative agreement with the TDSE result, except for very low energies. 
This shows that the shift of the interference maxima and minima should be en- 
tirely attributed to the effect of the Coulomb field. As we found in calculations, 
the most significant contribution to the shift is coming from the real part of (|13p 
induced by the Coulomb interaction in the Lagrangian. Since the action Wq is 
logarithmically divergent at the upper integration limit one should calculate it 
along an exact trajectory r±(p,t) where the asymptotic momentum is p. Only 
in this case the unphysical divergent parts cancel each other in the phase differ- 
ence of the interfering trajectories. For the calculation of the spectrum along the 
polarization direction, as it was done for Fig. 3, the trajectories i"o(p,t) cannot 
be taken as a zero approximation in the iteration procedure since half of them go 
exactly through the origin. According to Subsection 3.2, the proper trajectories 
fo should carry a nonzero transverse momentum. Therefore, trajectories which 
enter the iteration procedure must also have some nonzero transverse momen- 
tum. In this case the iteration procedure does not require a smooth turn-on of 
the nuclear charge (as it was done for elliptical polarization). The procedure 
yields rapidly converged trajectories, and the renormalized transverse momen- 
tum p y does not depend on the particular (non-vanishing) seed value for the 
transverse momentum. 

4.3 Other possible applications 

There are at least two other straightforward applications of the theory. 
1. With the development of intense sources of coherent XUV radiation based 
on free electron lasers delivering photon energies fouj ~ lOeV and higher [50j . 
the experimental study of strong field ionization of atoms and positive ions as 
well as other related processes in the intermediate (7 ~ 1) and multiphoton 
7 ^> 1 regimes became possible [51, 52, 153] . In particular, ionic charge states 
up to Xe 21+ were detected from a xenon gaseous target subject to 13.4 nm 
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Figure 3: (Color online) Photoelectron spectrum along the polarization direction 
calculated from the standard SFA (green), CCSFA (red) and the ab initio TDSE 
solution (black) for ionization from the ground state of hydrogen by a linearly 
polarized field with £ = 0.084 (2.5 • 10 14 W/cm 2 ) and uj = 0.056 (1.55eV). The 
TDSE solution is performed for an 8-cycle pulse (|26|) while the SFA and CCSFA 
amplitudes were calculated for the field IT8|) of constant envelope. The spectra 
are shifted vertically for visual convenience. 
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wavelength pulses of intensities about 10 16 W/cm 2 [SI]. Under such conditions 
a typical value of the Keldysh parameter is 7 = 10-^50. Even the single electron 
ionization dynamics is still poorly understood in this frequency domain. The 
method introduced in this paper can be applied to the evaluation of the Coulomb 
correction to the total ionization rate in the multiphoton regime, 7 > 1, i.e., 
for the above described conditions. This task was partially considered in [55] 
where it was shown that this correction may be even larger in absolute value 
than in the tunneling limit (where it is also substantial |24J. 
2. The method can also be used to calculate the probability of excitation into 
a Rydberg state, instead of ionization. Indeed, in short pulses the electron 
liberated from the atom and propagating towards a detector may have a negative 
total energy when the pulse is off. In this case such an electron will be captured 
in Rydberg states. The effect of anomalously high Rydberg population was seen 
in ab initio simulations and recently found in an experiment |56j . In order to 
obtain a distribution over the Rydberg levels one should apply the algorithm 
described above assuming a given negative energy and an angular momentum 
corresponding to a Rydberg state, instead of a final photoelectron momentum 
p corresponding to ionization. 

5 Conclusions 

The Coulomb corrections we introduced here are based on trajectories and are 
therefore gauge invariant. As a result the CCSFA amplitude is entirely gauge 
invariant for ionization from a s-ground state. For excited states, there is some 
gauge dependence appearing via the angular part of the bound state wave func- 
tion ^0 ( r ) m a short-range well we use to calculate the pre-exponential in ([24]) . 
This dependence is, however, much weaker than the one which appears in the 
standard SFA if an atomic wave function is being used in |T]). 

The application of our method to the strong field ionization of molecules 
seems to be the most natural generalization of the theory. Besides, per con- 
struction, CCSFA is applicable to any system where a long-range single particle 
or self-consistent force is important. Among such systems small metal clusters 
and fullerenes can be mentioned. Moreover, since in these systems the Coulomb 
interaction is not singular, evaluation of the CCSFA amplitude should be an 
even simpler task than it is for atoms and molecules. Finally, the method can 
be applied not only to ionization but also to other strong field processes involv- 
ing charged particles. 
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